Structure of the icosahedral Ti-Zr-Ni quasicrystal 
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The atomic structure of the icosahedral Ti-Zr-Ni quasicrystal is determined by invoking similarities 
to periodic crystalline phases, diffraction data and the results from ab initio calculations. The 
structure is modeled by decorations of the canonical cell tiling geometry. The initial decoration 
model is based on the structure of the Frank-Kasper phase W-TiZrNi, the 1/1 approximant structure 
of the quasicrystal. The decoration model is optimized using a new method of structural analysis 
combining a least-squares refinement of diffraction data with results from ab initio calculations. 
The resulting structural model of icosahedral Ti-Zr-Ni is interpreted as a simple decoration rule and 
structural details are discussed. 
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I. INTRODUCTION 

Quasicrystals are solids that combine long-range trans- 
lational order with rotational symmetries that forbid pe- 
riodic translational symmetry. To understand the phys- 
ical properties of quasicrystals, and in particular, the 
conditions under which nature prefers quasiperiodic to 
periodic order, it is necessary to determine their atomic 
structure. The aperiodic nature of quasicrystals requires 
new methods for structural analysis. In spite of signifi- 
cant progress in recent years, the complete atomic struct, 
ture of quasicrystals still remains an unsettled question.El 

The purpose of this paper is to determine the atomic 
structure of the thermodynamically stable icosahedral 
TiZrNi quasicrystal (i-TiZrNi). Although this material 
is the best-ordered among the known Ti-based quasicrys- 
tals, only polycrystalline samples with small grain sizes 
are available, which limits the information available from 
diffraction experiments. Therefore this paper develops a 
new method for structural analysis, based on the real- 
space structure description, whereby (several kinds of) 
tiles fill space and each kind of tile receives an identical 
"decoration" by atomic species, just as unit cells do in 
ordinary crystallography. 

The-.known quasicrystal structures fall into two main 
types.cl In a Frank-Kasper type structure, neighbor- 
ing atoms always form tetrahedra, and the constituent 
species divide into "small" atoms (with icosahedral coor- 
dination) or "large" atoms (with coordination number up 
to 16). In an Al-transition metal (Al-TM) structure, the 
transition metals are commonly surrounded by an icosa- 
hedron of Al, and Al atoms often form octahedra. In ei- 
ther type, icosahedral quasicrystals (and closely similar 
crystalline phases called "approximants" ) may be built 
from large (~50 atom) clusters with icosahedral symme- 
try: the "Bergman cluster"Ju in the Frank-Kasper case, or 
the "Mackay icosahedron"! in the Al-TM case. The Ti- 
based class of quasicrystals appears to fall on the border- 



line between the two main typesi; thus, though i-TiMn 
and i-TiCr fall in the Al-TM class, our subject i-TiZrNi 
is of the Frank-Kasper type. In developing its structural 
model, we will take advantage of structural relationships 
to both the i-AlMnSi and the Frank-Kasper quasicrys- 
tals. 



A. Methods of Structural Investigation 

Many thermodynamically stable quasicrystals such as 
i-AlCuFeO, i-AlPdMnfl, and d-AlMCoQ within the Al- 
metal class, i-AlCuLi and i-MgZnRE (RE=rare earth)£ 
in the Frank-Kasper class, as well as the new binary qua- 
sicrystal i-CdYbl!2l, can be grown as single grains up to 
several millimeter size, which allows the measurement 
of single crystal diffraction data. (Of those mentioned, 
all except «-AlCuLi also exhibit long-range order, i.e., 
instrumentally. sharp Bragg peaks.) Using a variety of 
proceduresfilOrB, atomic structures have been tefpcd. for 
many oL-thfise quasicrystals, pich-as i-AlCuLOolia, i- 
AlCuROiMS and i-AlPdMnBBB. 

All these structures were parametrized as cuts by a 
three-dimensional hyperpi^pcy-through a density in a six- 
dimensional hyperspace.c3'CJ'cJ The difficulty with the 
six-dimensional approach is that only an average struc- 
ture can be determined and disorder in the quasicrystal 
can not be treated appropriately. This, combined with 
experimental limitations such as truncation effects, usu- 
ally produces a small number of spurious atoms with un- 
physical coordination numbers or small interatomic sep- 
arations in structural models. A purely crystallographic 
approach alone might, therefore, not be able to fix all the 
positions and chemical identities of the atoms in qua- 
sicrystalline structures. It is furthermore infeasible for 
quasicrystalline alloys with small grain sizes for which 
only powder diffraction data are available. 

Additional experimental and theoretical data, such as 
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structural information from related approximant phases 
and results from total energy calculations, can be helpful 
in the development of realistic structure models for qua- 
sicrystals. To incorppcate this information a real space 
approach is desirablerir 6 ! A common real space approach 
to the structure of quasicrystals is the tiling decoration. 
One splits the structural analysis into two parts: (1) the 
problem of the determination of the correct tiling ge- 
ometry, and (2) the atomic decoration of the tiles. An 
advantage of the decoration viewpoint is different calcu- 
lations (ab initio LDA, pair-potentials, diffraction simu- 
lation) may be performed using periodic tilings (" approx- 
imants") with larger or smaller unit cells, depending on 
the calculation's computational burden, and the results 
can be compared or combined in a systematic fashion. 

The first tiling-decoration approaches wereJpased on 
similarities to periodic crystalline structures. In later 
work energetic considerations were taken into account. 
Examples of such investigations are the work by Krajci 
and Hafner who investigated the structure of the i- 
AlMgZn quasicrystal using a decoration of the three- 
dimensional Penrose tiling and pair potentials based on 
the pseudopotential theory.EZl Windisch et al. studied 
the atomic and electronic structure and the lattice dy- 
namics of i-AlCuLi using molecular dynamics simula- 
tions and density-functional calculations of decorated 
Penrose and canonical cell tilings P^P" The structure of 
the (metastable) i-AIMn quasicrystal has been studied by 
Mihalkovic et al. usipg canonical cell tiling decorations 
and pairpotentials.E§E3 Recently Mihalkovic et al. inves- 
tigated the structure of decagonal AlCoNi using Monte- 
Carlo calculations with tile updateslill 

All these investigations are based on the energetics of 
structural models and do not incorporate diffraction data 
into the modeling process. To derive a reliable structural 
model of a quasicrystal, a combination of the crystallo- 
graphic approach based on diffraction data with informa- 
tion from total energy calculations and relaxation stud- 
ies is highly desirable. In this paper we present the first 
study combining a least-squares refinement of X-ray and 
neutron diffraction data with structural information ob- 
tained from ab initio relaxations for approximant struc- 
tures. 



B. The Ti-Zr-Ni quasicrystal 

Few studies have been performed for quasicrystalline 
alloys with grain sizes so small that only powder diffrac- 
tion data are available. Earljtjijoxk was devoted to the 
metastable system i-AlMnSi.EjOEfl Another class of such 
quasicrystals are the Ti-based icosahedral alloys. As 
early as 1985, Zhang et al.E3 discovered a metastable 
icosahedral phase in Ti-Ni-V. In the following years, 
further icosahedral phases in Ti-based alloys containing 
other 3d transitk^mfttak^ vLch. as V, Cr, Mn, Fe, Co and 
Ni were found.^OH'^c3cllc3 However, no stable qua- 
sicrystalline phase of this class has been discovered yet. 



Investigations of projection-based models of the struc- 
ture of the Ti-Mn quasicrystal by Phillips et al. using 
tight-bindingnxecursion and pair-potential based energy 
calculationsocJ showed that the Mackay clusters form 
the backbone of the structure. Furthermore small ad- 
ditions of Si and O, which are cruciaL£ar the formation 
of the Ti-Cr and Ti-Mn quasi crystaJocM, were found to 
energetically stabilize the structureo 

In this work the experimentally well-studied icosahe- 
dral Ti-Zr-Ni quasicrystal is considered. As early as 1990, 
Molokanov and ChebotnikovCJ discovered a metastable 
icosahedral phase in Ti-Zr-Ni and in 1997 a thermo- 
dynamically stable icosaiedral Ti-Zr-Ni phase was dis- 
covered by Kelton et al.cS The i-TiZrNi quasicrystal is 
formed by either rapid quenching or solid state trans- 
formation at temperatures of 500°C to 600°C, generally 
leading to a fine microstructure of quasicrystal and crys- 
tal phases. The icosahedral phase is rather disordered 
with a coherence length of about 350 A, resulting in a 
small number of measurable structure factors. Therefore, 
an atomic structure determination by diffraction experi- 
ments alone is impossible. Consequently little is known 
about the atomic structure of i-TiZrNi. 

There is hope of obtaining information on the atomic 
arrangement in i-TiZrNi by studying the structure of 
related periodic "approximants" , meaning crystal struc- 
tures whose unit cells are fragments of some quasicrystal. 
There are several competing periodic phases known in 
the Ti-Zr-Ni phase diagram. Besides the binary phases 
a-TiZr, <5-Ti2Nio and /-Zi^NiE^, there are two ternary 
phases; a cubic Frank-Kasper type structure, (VT^-phase) 
and a hexagonal Laves type phase A-TiZrNi.El The struc- 
tures of the VF-TiZrNi phase, which can be interpreted 
as a Fibonacci 1/1 approximant of i-TiZrNi, and of the 
A-TiZrNi phases were determined by a Rietveld refine- 
ment using X-ray diffraction and neutron scattering data. 
Struc tura l d etails were confirmed by ab initio calcula- 
tions .e3Ejo The Wf-TiZrNi phase is closely related to the 
quasicrystal phaseH£3 They both form at similar com- 
positions and temperature ranges and grow with a coher- 
ent crystallographic orientation. EXAFSrStudies indicate 
a similar local structure in both phases.LB The presence 
of a reversible phase transformation between the icosa- 
hedral phase and the VF-phase at about 600°C indicates 
that z-TiZrNi is a low-tepajperature phase, in contrast to 
most other quasicrystals£3 This opens the exciting possi- 
bility that an (aperiodic) i-TiZrNi quasicrystal, or a very 
similar large-cell crystal, could be a ground state phase. 



C. Contents of this paper 

The goal of our work is to model the ideal structure 
of the i-TiZrNi quasicrystal. The energetic stability and 
electronic structure of approximants to the «-TiZrNi qua- 
sicrystals, as well as the competing phases, are the topics 
of a subsequent paper. In this paper we present a struc- 
tural analysis of z-TiZrNi based on a new method, com- 
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bining powder diffraction data and information from ab 
initio calculations in form of a constrained least-squares 
analysis. In Section II the initial decoration model is de- 
rived. We base our structure on the canon ical cell tiling 
which is 
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geometry, which is introduced in Section 
starting point an initial atomic decoration of the canon- 
ical cells is determined from the at omic structure of the 
1/1 approximant IF-TiZrNi (Sec. [ilCf) . The following 
Section III on the optimization of the decoration model 
is the core of the paper. The constrained least-squares 
method of refinement is introduced. This approach in- 
corporates structural information from ab initio relax- 
ations as constraints into a least squares refinement of 
the diffraction data. The resulting structure can be un- 
derstood as a simple decoration model. This and further 
details of the structural model, which are important for 
comparisons to experimental data such as EXAFS mea- 
surements, are discussed in Section IV. 



II. FIRST DECORATION MODEL 

In this section we derive the initial decoration model. 
The underlying tiling geometry of our model is given by 
the canonical cell tilings. The canonical cell tilings pro- 
vide a tiling model with a high packing fraction of icosa- 
hedral clusters. The location of the atomic sites dec- 
orating the canonical cells in our model are based on 
structural features of the VF-phase. 



TABLK-iI: Crystallographic description of the phase W- 
TiZrNiLJ showing the chemical occupation of the lattice sites. 
The positions of the sites are given in units of the lattice con- 
stant. 
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Ti(1.0) 


8 


Pi 


0.412 





0.5 


Zr(0.5), Ti(0.5) 



The atomic structure of the IF-phase is similar to the 
structure of tha Frank-Kasper phases Al4oMg4oZn2otl and 
Al 56 Li32Cui 2 Ea. The IF-TiZrNi structure differs from 
these in (i) the occupation of the central site of the 
Bergman cluster, which is empty in the cases of Al-Mg- 
Zn and Al-Li-Cu, and (ii) the chemical disorder, which 
occurs between the smaller components Al and Zn or 
Cu in Al-Mg-Zn and Al-Li-Cu respectively and between 
the larger components Ti and Zr in the case of Ti-Zr-Ni. 
Furthermore, the analysis for i-TiZrNi showed chemical 
ordering on the sites of the Bergman cluster. 



B. Canonical cell tiling 



A. Structure of the W-phase 

The VF-TiZrNi phase is a periodic crystalline phase 
closely related to the icosahedral quasicrystal. Its com- 
position TisoZrssNils is close to the composition of the 
quasicrystal Ti4i.5Zr4i.5Nii7 and it grows with a coher- 
ent crystallographic orientation with respect to the qua- 
sicrystalline phase. The structure of the VF-phase was 
determined using X-ray diffraction and neutron .scatter- 
ing data and confirmed by ab initio calculations E3E3 

It consists of a body centered cubic packing of icosa- 
hedral Bergman clusters connected by a few "glue" sites 
(see Table [[]). The Bergman cluster as shown in Fig- 
ure [l] is composed of a central Ni atom surrounded by a 
small icosahedron of 12 Ti atoms. A second icosahedron 
is formed by placing 12 Ni atoms directly outside the 
vertices of the small icosahedron. Finally 20 Zr atoms 
are placed on the faces of the large icosahedron. The re- 
maining 36 "glue" sites consist of 24 Ti atoms, forming 
a hexagon that separates the Bergman clusters along the 
threefold (111) directions, and 12 Ti and Zr atoms, form- 
ing a rhombus that divides the Bergman clusters along 
the twofold (100) directions. The structural analysis and 
ab initio calculations revealed that chemical disorder in 
the VF-phase is abundant only between the chemically 
similar Ti and Zr atoms on the "glue" sites but no dis- 
order is found on the cluster sites.Ea 



In modeling the quasicrystal structure the canonical- 
cell tiling (CCT) model is used.EJ The packing rules 
permit periodic, random and in principle quasiperiodic 
tilings, although ..no rule for a quasiperiodic tiling has 
been found yet.t3E3 A large number of periodic approx- 
imants structures, ranging from a 1/1 up to a 31/24 ap- 
proximant, are known for this tiling. In order to use the 
CCT's for a structural description of quasicrystals and 
for a refinement of decoration models to experimental 
diffraction data, predictions for the structure factors of 
the icosahedral structure are needed (see App. |X]and|^). 

Canonical cell tilings consist of four kinds of disjoint 
cells denoted A, B, C and D. Their nodes are joined 
by two types of linkages, b and c bonds, which run in 
twofold and threefold symmetry directions respectively 
(see Fig. ^). There are three different kind of cell faces. 
The X face has the shape of an isosceles triangle formed 
by one b bond and two c bonds, the Y face with a shape 
of an equilateral triangle is formed by three c bonds and 
finally the rectangular Z face is formed by two b bonds 
and two c bonds. 

The simplest quasiperiodic tiling with icosahedral sym- 
metry is the Penrose tiling. Its tiles are the prolate rhom- 
bohedron (PR) and oblate rhombohedron (OR); their 
edges have the length a q and run into fivefold symme- 
try directions. The CCT can uniquely be decomposed 
into PR's, OR's and rhombic dodecahedra consisting of 
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two PR and two Oi?.E3E3 A rhombic dodecahedron is 
assigned to each b bond and a PR to each c bond. They 
are called RD and PR C respectively (see Fig. ||). They 
both connect nodes of the CCT. 

Each triangular Y face is pierced by another kind of 
PR, which we will call PRy. It is useful to divide these 
PRy into two subclasses, called PRy(C) and PRy(D) 
according to whether it lies mostly in a C or D cell re- 
spectively. Each D cell contains three more Pi?'s, called 
PRd- Finally, each Z face contains an OR. This de- 
composition of the CCT facilitates the description of the 
decoration model in the following section. 



C. Decoration description 

We now derive the initial decoration sites for the 
canonical cell tiles from the structure of the W-phase. As 
our starting point we require that our decoration model 
reproduces the features of the VF-TiZrNi structure wher- 
ever possible. 

The structure of the VF-phase corresponds to a Fi- 
bonacci 1/1 approximant. In the canonical cell tiling 
description the 1/1 approximant is given by an Aq tiling. 
The decoration of the A-cell is thus completely deter- 
mined by the VF-phase. The nodes of the Aq tiling cor- 
respond to the center of the Bergman cluster in the W- 
phase. Thus every node of the CCT will be occupied by a 
two-shell Bergman cluster and the atomic configuration 
surrounding the twofold (b) and threefold (c) linkages will 
be the same as in the W-phase. This already accounts for 
over 80% of the atoms in decoration of large tilings. We 
then extend our decoration model to the interior of the 
larger B, C and D cells using the structure of W-TiZrNi 
as a guide. 

We use the decomposition of the larger canonical cell 
tiles into prolate and oblate rhombohedron as well as 
the rhombic dodecahedron to describe their decoration. 
Since the Aq tiling can be decomposed into Pi? c 's and 
P-D's, their atomic decoration is determined by the 
atomic structure of the 1/1 approximant (see Fig. 
Simply to limit the number of free parameters, it is as- 
sumed that the decoration sites of the PRy(C), PRy(D) 
and PRd are the same as for the PR C . As for the OR, 
its face decoration is forced by the need to adjoin with 
the other tiles and it has no space for additional interior 
atoms (see Fig. ||). It is noted here that the site list is the 
same as for the Henley-Elscr decoration of Al-Mg-Zn.H 

After defining the general sites of the decoration model, 
the nomenclature for the different sites, necessary for 
a description of the chemical decoration, is introduced. 
From the VF-phase we derived four general types of dec- 
oration sites for the Amman tiles. These are the vertices 
and the mid-edge sites of the Amman tiles, the two diag- 
onal sites of the PR, and the sites inside of the RD. We 
denote these four kinds of sites by the letters a, (3, 7 and 
8 respectively (see Fig. ^). Any further differentiation of 
the sites is indicated by a subscript. The physical sig- 



nificance is that identically labeled sites have the same 
coordinates and occupations in the decoration model. 

The center and the second shell icosahedron of the 
Bergman cluster decorate vertices of the Amman tiles; 
these are denoted ot\ and ai respectively. The sites of 
the first shell of the Bergman cluster are equivalent to 
edge centers of the Amman tiles and are denoted by 
The PR C has six mid-edge sites on the plane bisecting 
its long axis. The two atomic sites on the long axis of the 
PR C which, furthermore, also decorate the faces of the 
large icosahedron of the Bergman cluster, are denoted by 
71 . The same 71 site also appears in the interior of the 
RD. There are four additional sites in the interior of the 
RD on the plane bisecting its long axis. These four sites 
form a rhombus. The sites on its short axis are denoted 
by Si and the sites on its long axis by 62- All of the dec- 
oration sites discussed up to now were directly derived 
from the structure of the W-phase. 

For the B, C and D cell we assume that the local struc- 
ture is similar, with the decoration of the nodes and the 
b and c bonds the same as derived from the M^-phase. 
There are additional sites which need to be specified for 
the B, C and D cell. These are the sites connected to 
the PRy(C), the PRy(D), the PR D and the OR. The 
acute vertex of the PRy closest to the Y face, the three 
mid-edges sites adjacent to this vertex and the adjacent 
site on the threefold axis are not already determined by 
adjacent tiles. These sites are called 0:3, 0s and 72 re- 
spectively (see Fig. |5](a)). For the PRy(D), a number of 
sites inside the D cell must be specified. Starting from 
the acute vertex away from the Y face, the three mid- 
edge sites /?5, the three vertices 014, and the ring of six 
mid-edge sites need to be specified. Furthermore, the 
site along the diagonal that is adjacent to the acute ver- 
tex that is away from the Y face is denoted by 73 (see 
Fig. ||(c)). For the OR, the mid-edge sites are labeled 
/?4 (see Fig. ||(d)). Along the diagonal of the PRd, the 
site away from the node of the D cell is called 74. The 
six mid-edge sites on the plane bisecting the PRd are 
denoted (3 7 (see Fig. |(b)). 



III. STRUCTURAL REFINEMENT 

In this section the decoration model for i-TiZrNi is 
optimized using a constrained least-squares analysis of 
diffraction data incorporating results from ab initio cal- 
culations. In this approach, the optimization of the struc- 
tural model is split into two steps. In the first step, the 
chemical decoration of the atomic sites is refined by X- 
ray and neutron diffraction data while the site positions 
are kept fixed. The atomic positions are then optimized 
in a second step by ab initio relaxations of the result- 
ing atomic decoration for periodic tilings. These relaxed 
atomic positions are then used as a constraint to the re- 
finement of the diffraction data. 
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A. Diffraction data 

The sample preparation and the data acquisition was 
performed by E. H. Majzoub; details can be found in 
References 
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The composition of the sample is 
(Ti4i. 5 Zr4i. 5 Nii7798Pb2. The Pb was used, since it was 
assumed to accumulate in the grain boundaries due to 
its low solubility in the Ti-Zr-Ni phases. The regions 
to which the Pb segregates are liquid during the anneal- 
ing, allowing for faster diffusion of the transition metals, 
thus enhancing the formation of larger quasicrystallinc 
grains with longer coherence length. The samples were 
arc-melted from the pure materials under an argon at- 
mosphere and then annealed at 600° C for 7 days under 
vacuum. The X-ray diffraction pattern was measured 
using Cu-K Q radiation. Neutron powder diffraction data 
for the sample was obtained at the University of Missouri 
Research Reactor in collaboration with W. B. Yelon. The 
neutron wavelength was 1.7675A. 

The uncorrected structure factors were obtained from 
the powder diffraction data by performing a least-squares 
fit of the data to a Gaussian peak profile, 



I{9) =I b + I 



1 



2ttct 



exp 



2<j 2 



where the peak is centered at 6q and has a width of 
a. The integrated intensity is given by 1$. Any peak 
broadening related to the diffractometer or the coher- 
ence length of the sample is absorbed in the width, a, 
of the Gaussian peak profile. The data were fitted in a 
small interval of AO f» 1° ... 2° around the center of each 
peak. The dependence of the results on the choice of the 
size of this interval was accounted for by the estimated 
errors, which reflect the range over which the fit param- 
eters varied when the size of the interval over which the 
fit was performed was slightly changed. The background 
was fitted using a single parameter, lb. Out of the 31 
peaks for the neutron diffraction data, 24 peaks were in- 
dexed to the icosahedral phase and 7 were indexed to the 
(5-Ti2Ni secondary phase. For the X-ray data, 18 struc- 
ture factors were obtained. Three of the peaks belonged 
to the <5-Ti2Ni secondary phase and 15 to the icosahedral 
quasicrystal. 

Overall, the number of measured structure factors is 
rather small. Furthermore, since only powder diffrac- 
tion data is available for the i-TiZrNi quasicrystal, a few 
of the measured structure factors were superpositions of 
crystallographically inequivalent reflections, further ex- 
acerbating a structural analysis. In the following two 
sections the initial decoration model for canonical cell 
tilings is refined to the measured structure factors con- 
straining the atomic positions by ab initio relaxations. 



B. Calculating the structure factors 

A least-squares analysis for the chemical decora- 
tion of the sites of the structure model to the X-ray 



diffraction and neutron scattering data from icosahedral 
Ti41.5Zr41.5M17 was performed. As a good approxima- 
tion of the quasiperiodic tiling, a cubic 8/5 CCT approx- 
imant with 576 nodes forming 1320 A, 576 BC and 136 
D cells containing 52392 atoms per unit cell and with a 
lattice constant of 98A was used.o For the 8/5 approxi- 
mant, the distribution of the complementary coordinates 
of the nodes has a small variance. It is thus reasonable 
to assume that this tiling is a good approximant to the 
hypothetical quasiperiodic canonical cell tiling. The dec- 
orations of the 8/5 canonical cell tiling were constructed 
using the Canonical Cell, Tiling Package - "There are the 
atoms" by MihalkovicS 

The icosahedral structure factor, i*i co (q), is approxi- 
mated by averaging over the amplitudes, -FkpplQa), of all 
the Bragg peaks in the approximant that map onto a sin- 
gle orbit of icosahedral peaks given by the wave vector, 

q, 



F i<=o(q) = 1 ^M ap p(q a )F ap p(q a ). 
-"'-'icol.qj a 



(2) 



Each corresponding orbit of reflections in the approx- 
imant is represented by a wave vector, q a , and has the 
(1) multiplicity M app . M; co is the multiplicity of the icosahe- 
dral orbit. The average is taken over amplitudes instead 
of intensities requiring the correct choice of phases of the 
structure factors. It is therefore necessary to center the 
approximant structure in six dimensional space. Details 
of this and a discussion of the differences between av- 
eraging over amplitudes and intensities can be found in 
Appendices |X] and 

The X-ray and neutron diffraction data were taken at 
room temperature. To account for thermal motion which 
decreases the structure factor, Fi CO (q), for large momen- 
tum, q, an isotropic Debye- Waller factor is used, 



/dw(?) = ex P 



D5 

.^igo 2 

8ir 2 q 



(3) 



Due to the small number of measured structure fac- 
tors, only a single thermal parameter, B? so with s S 
{Ti, Zr, Ni}, for each atomic species was used. The X-ray 
form factors, f s (q), and the neutron scattering lengths, 
f s , were taken from Reference Q . 

Altogether the icosahedral structure factors, Fj co , 
needed for the refinement are given by 



Fico(q) 



E^(o/Dw( g )/,( ? )E Mapp(qa) 



I,. 



M ico (q) ^ 



(4) 

The first sum is over all classes, I, of distinct sites of the 
decoration model. The second sum is over the three dif- 
ferent atomic species, s G {Ti, Zr,M}, decorating these 
classes of sites with concentrations, c s (l). The third sum- 
mation is over the orbits of the approximant structure 
represented by the wave vector, q a , which correspond to 
the icosahedral wave vector, q. Finally the last summa- 
tion is taken over the individual atomic site, Vj , belonging 
to the set of sites, Ci, of the site class, I. 
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C. Constrained least-squares refinement 

Before we can compare the experimental peak intensi- 
ties to the calculated structure factors, several correction 
factors have to be taken into account. For both, the 
X-ray and the neutron diffraction data, the intensities 
must be corrected by the Lorentz factor, sec cosec 2 0. 
Furthermore for the X-rays, the intensities are modified 
by the polarization, P, of the beam as described by the 
Thomson formula 



1 + Pcos 2 (20) 



P 



(5) 



Since the X-ray beam was unpolarized, P = 1 is used. 
Overall the estimated intensity of a reflection is given by 



-fico(g) = 



1 + Pcos 2 (20) 



(f + P)(cos0sin 2 0) 
for X-ray diffraction and by 

^ico(g) = sec0 cosec 2 |P ico (g)| 2 



l^ico(«)r (6) 



(7) 



for neutron scattering. 

Since the positions of the atoms are kept fixed in the 
refinement, the last two summations in Equation (|]) over 
the atomic sites of each class with identical atomic deco- 
ration and over the orbits of wave vectors for the approx- 
imant are performed beforehand for all necessary values 
of qi co . The total structure factor for each q was then 
obtained by summing over these partial structure fac- 
tors multiplied by the average form factor or scattering 
length, f s (q), and the average Debye- Waller factor of the 
atoms decorating the site. The corrected intensity, 7i CO , is 
then obtained from the calculated structure factor, F- lco , 
using Equations || or [7] respectively. 

From the calculated intensities, the residuals R and R w 
as well as the reduced x 2 are obtained, providing a stan- 
dard measure of how well the calculated data matches 
the observed one. The residuals are defined as 



R 



E 9 VW?) - \Aco(?) 



(8) 



and 



Rqn — 



\ 



Eg W (?) J exp(9) 



(9) 



where the sums are over all observed and calculated 
intensities, I cxp and ii co respectively. The weights w 
are given by the uncertainty of the measured intensities 
where w — 1/cf. The "goodness" of the fit is given by 
the reduced x 2 , 



Eg w(q) (yr^(q) - \Aco(<?)) 

N — N n 



(10) 



where N is the total number of reflections and N p is the 
number of estimated parameters. 

Since the number of measured reflections is rather 
small, additional measured quantities such as the compo- 
sition of the structure can be incorporated as constraints 
to the refinement and to limit the parameter space. The 
composition of the refined structure, a, was constrained 
by a penalty function, 



ie{Ti,Zr,Ni} 



c°)> 2 



(11) 



In order to ensure that the composition of the refined 
structure, c,, does not deviate by more than a few percent 
from the nominal experimental value, c°, a value of a c = 
0.02 was used. The \ 2 was then minimized by varying the 
chemical decoration of the sites using Powell's conjugate 
direction method. E3 

The chemical decoration of the ideal atomic sites is re- 
fined with respect to the structure factors of 38 reflections 
for i-TiZrNi, 15 X-ray and 23 neutron diffraction struc- 
ture factors. The maximum number of different sites of 
the decoration model is 18. The decoration of the sites 
in the W-phase structure was used for the initial config- 
uration, i.e. a sites were decorated with Ni, /3 sites with 
Ti and 7 and S sites with Zr. The chemical decoration 
on each of these site classes was allowed to vary between 
two components, either Ti and Zr or Ti and Ni. About 
a dozen different starting configurations were used. The 
majority of the optimization runs resulted in the same 
final concentrations and gave the lowest observed x 2 ■ 
However, as common in optimization problems, we can- 
not guarantee to have found the global minimum. Dur- 
ing the refinement classes of sites which showed similar 
chemical occupancies were combined into a single class. 
The final number of independent classes of sites was six. 
This, and the three thermal parameters yield a total of 
nine estimated parameters. 

The refinement using the initial decoration sites 
yielded for the residuals, R, of the structure factors 
11.5% for X-ray, 6.1% for neutron and 6.9% combined. 
The weighted residuals, R w , are 15.7% for X-ray, 4.6% 
for neutron and 5.7% combined. The reduced x 2 is 
1.8 indicating crystallographic significance for the fit. 
The thermal parameters for Ti and Zr, B^ a — 1.1 A 2 
and Bg = 1.2 A 2 , are_comparable to values found in 
the <5-(Ti,Zr) 2 Ni alloys.B The thermal parameter for Ni, 
-E?iso = 2. 3 A 2 , is about a factor two larger than the value 
found in <5-(Ti,Zr)2Ni. This could reflect static disorder 
on some of the Ni sites. 



D. Ab initio relaxations 

So far we have only refined the chemical decoration of 
the ideal sites of the canonical cell tiling. In the next 
step we perform ab initio relaxations for the resulting 
structure to investigate the site positions. The ab initio 
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calculations were performed with VASPl 



which is 



a density functional code using a plane-w, a> r i e .basis and 



Bfi 



ultrasoft Vanderbilt-type pseudopotentialsLJ'ILU. Atomic- 
level forces are calculated and relaxations with a con- 
jugate gradient method are performed. Ab initio calcu- 
lations are computationally very expensive. Thus, only 
calculations for small systems with up to a few hundred 
atoms are feasible on current computers. Therefore, ab 
initio relaxations only for the small periodic tilings of the 
canonical cells were performed. Tie tilings used were the 
A 6 , the B 2 C 2 and the D 2 tiling E3 

As a starting configuration, the results from the analy- 
sis of the diffraction data were used. The chemical order 
was idealized by placing only the majority species on each 
site. The calculations were performed using the general- 
ized gradient approximation by Perdew and Wang£3, and 
a plane- wave kinetic-energy cutoff of E cut = 302.0 eV was 
chosen to ensure convergence. The pseudopotentials for 
Ti and Zr describe, besides the usual outer shell states, 
also the 3p and 4p states, respectively, as valence states. 
Due to the size of the quasicrystal-like tiling structures a 
check of the convergence with the number of k-points 
was not always possible due to computational limita- 
tion. Therefore the k-point mesh was chosen as large 
as possible, i.e. 3 x 3 x 3 for the A 6 tiling, 2x2x2 
for the B2C2 tiling and 1 x 1 x 1 for the D 2 tiling, 
corresponding to 4, 2 and 1 k-points in the irreducible 
part of the Brillouin zone respectively. The positions of 
the atoms, as well as the shape and volume of the unit 
cells, were relaxed until the total electronic energy change 
was smaller than 1 meV. This corresponds to atomic-level 
forces F max < 0.02 eV/A. 



TABLE II: Maximum and average displacements in [A] of 
the atoms from their ideal positions for the different periodic 
canonical cell tilings after relaxation. 



Tiling Atoms Displacement 

average (rms) max 



A 6 
B2C2 
D 2 



81 
91 
123 



0.10 
0.11 
0.13 



0.19 
0.22 
0.25 



The atomic positions after the ab initio relaxation were 
remarkably close for all three tilings to the ideal atomic 
positions, such as the corners, mid-edges etc. of the dec- 
orated tiles, with average displacements smaller than 
0.13 A (see Table 0). 

The resulting atomic positions from the relaxation cal- 
culations were then symmetrized. This involves mapping 
the relaxed positions back onto their canonical orbit and 
calculating the, average canonical coordinates of equiv- 
alent atoms.E3lEl The deviations of the relaxed atomic 
positions from their canonical coordinates are in general 
rather small with an average deviation of 0.07 A. The 
largest deviation occurs for the ft site, where the relaxed 
positions deviates by 0.19 A in the D 2 and by 0.05 A in 



TABLE III: Results of the refinement of the chemical decora- 
tion for the quasicrystal using the relaxed atomic positions. 



Site 


Description Decoration 


Class 




Cluster 


Center 


Ni (1.0) 


1 
1 


ft 




Small icosahedron 


Ti (1.0) 




A 


0.1 




Large icosahedron 


Ni (1.0) 


1 
± 


71 




r aceo 


7r (1 n\ 


3 


Si 


6-bond 


RD Close pair 


Zr (1.0) 


3 


62 




RD Distant pair 


Zr (1.0) 


3 


ft 


c-bond 


PR C Mid edge 


' Ti (0.8) 

7,r (ft 9} 


4 


Ot3 


Y-face 


PRy(C) Vertex 


Ni (1.0) 


1 


ft 




PRy{C) Mid-edge 


Ti (1.0) 


2 


72 




PRy(C) Upper diagonal 


' Ti (0.6) 
Zr (0 4\ 


5 


ft 


Z-face 


OR Mid-edges 


' Ni (0.7) 
^ Ti (0.3) 


6 


73 
Q4 


D-cell 


PRy(D) Lower diagonal 
PRy(D) Vertices 


Zr (1.0) 
Ni (1.0) 


3 
1 


ft 
74 
ft 
ft 




PRy(D) Mid-edges 
PRd Diagonal 
PRd Mid-edges 
PRd Mid-edges 


Ti (1.0) 
Ti (1.0) 
Ti (1.0) 
Ti (1.0) 


2 
2 
2 
2 



the B 2 C 2 cell. In fact, the deviations in both cells have 
almost opposite directions. 

The symmetrized positions are used as an input for 
fitting the site chemistry to the diffraction data. This 
process can, in principle, be repeated until convergence is 
achieved. In our case, since no significant changes in the 
decoration were observed, the optimization was already 
converged after the refinement of the relaxed atomic po- 
sitions. 

Using the relaxed positions, the residuals, R, of the 
structure factors decrease to 10.8% for X-ray, 5.1% for 
neutron and 5.9% combined. The weighted residuals, 
R w , are 15.0% for X-ray, 4.0% for neutron and 5.1% com- 
bined. The reduced x 2 is 1.4, indicating a good agree- 
ment between the calculated and measured structure fac- 
tors (see Fig. ||). Overall, the relaxed positions lead to 
a better agreement with the diffraction data than the 
unrelaxed positions. 

The resulting chemical decoration is given in Table III . 
The thermal parameter decreases for Ti to 0.3 A 2 , for Zr 
to 1.0 A 2 and it is unchanged for Ni. The chemical dec- 
orations on the sites change by less than 5%, with the 
exception of site 72 , where the probability of finding Ti 
increases by 10%, and the site 0:4, that is now only occu- 
pied by Ni. The stoichiometry changes by less than 0.5%. 
The composition of the refined structure Ti4iZr4oNiig 
agrees with the nominal composition Ti4i.5Zr4i.5Nii7 to 
within 2%. 

This is the first refinement that systematically com- 
bines diffraction data and structural energies. Since the 
powder diffraction data for i-TiZrNi are not available for 
sufficiently high momentum q to usefully constrain the 
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coordinates the atomic positions were fixed to the results 
from ab inito relaxations. If better data were available for 
this quasicrystal, positional information would be avail- 
able from the diffraction data as well. The relaxed atomic 
positions lead to an improved agreement between the cal- 
culated and measured structure factors. This illustrates 
that results from ab initio relaxations can successfully be 
incorporated as constraints into structural refinements. 



IV. STRUCTURAL FEATURES OF THE 
QUASICRYSTAL 

In this section we discuss some details of the resulting 
structural model for i-TiZrNi using the large 8/5 approx- 
imant and compare to the structure of the VF-phase. Re- 
sults for pair correlation functions, coordination numbers 
and local atomic environments, which are important for 
comparisons to exp erimental data such as EXAFS, are 
given in Tables |V]and |v| as well as Figure ^ 



TABLE IV: Number and average distance of nearest neigh- 
bors in parentheses in [A] in the VK-phase (a=14.3A). 



Center 


Ti 


Zr 


Ni 


Total 


Ti 


4.91 (2.90) 


5.12 (3.07) 


2.00 (2.66) 


12.03 (2.93) 


Zr 


7.97 (3.07) 


1.76 (3.12) 


2.75 (3.19) 


12.48 (3.10) 


Ni 


6.37 (2.66) 


5.63 (3.19) 


0.00 (-) 


12.00 (2.91) 



TABLE V: Number and average distance of nearest neigh- 
bors in parentheses in [A] in the 8/5 approximant model 
(a=98.0A). 



Center 


Ti 


Zr 


Ni 


Total 


Ti 


3.95 (2.88) 


5.85 (3.04) 


2.34 (2.67) 


12.14 (2.92) 


Zr 


5.91 (3.04) 


3.30 (3.02) 


2.91 (3.13) 


12.12 (3.06) 


Ni 


4.93 (2.67) 


6.05 (3.13) 


0.83 (2.70) 


11.81 (2.91) 



Figure |7j shows the partial pair-distribution functions 
for the 1/1 approximant phase, VF-TiZrNi, and the 8/5- 
approximant. Tables [fy| and |y| list the numbers and av- 
erage distances of nearest neighbors in both structures. 
The pair-correlation functions for the IF-phase and the 
quasicrystal resemble another, corresponding to the sim- 
ilar local atomic structures of both phases. While the W- 
phase has no Ni-Ni nearest neighbors at all, the number 
of Ni-Ni neighbor pairs in the quasicrystal is rather small. 
Chemical disorder was found only on a small number of 
sites, apparently localized to the glue region between the 
Bergman cluster, similar to the VF-TiZrNi approximant 
(see Tab. |). Results of recent investigations of the lo- 
cal atomic structure of the-W -phase and the quasicrystal 
by EXAFS measurementsE3 are consistent with the pair- 
distribution functions of the structural model developed 
here. 



The features of the decoration model of the quasicrys- 
tal can be understood by the energetics of the interac- 
tions between the constituents. Ti and Zr atoms are 
completely miscible and exhibit a zero heat of mixing. 
Ti and Ni as well as Zr and Ni, on the other hand, have 
a strongly negative heat of mixing, indicating strong at- 
tractive interactions between these pairs. Thus, it is en- 
ergetically favorable for all the Ni atoms to be surrounded 
by Ti or Zr. This explains why, in the structure of the 
1/1 approximant and the quasicrystal, hardly any Ni- 
Ni bonds occur. Since Zr is slightly larger than Ti it is 
no surprise that Zr occupies the more open sites of the 
structure given by the interior of the PR and RD. 

The refined structure of the quasicrystal is very sim- 
ilar to the structure of the 1/1 approximant phase, W- 
TiZrNi. The atomic sites for the structure model of the 
quasicrystal are based on the VF-phase structure. The re- 
fined decoration of these sites is the same in most cases. 
However, the Si and 82 sites along the b bond are occu- 
pied by Zr only in the quasicrystal, while in the IF-phase 
chemical disorder occurs between Ti and Zr. Overall, 
the amount of chemical disorder in the quasicrystal ap- 
pears to be of the same order or even slightly less than 
in the VF-phase. Considering the apparent disorder in 
the powder diffraction patterns, the lack of chemical dis- 
order on most sites of the structural model of the qua- 
sicrystal seems surprising at first and might be fortuitous. 
However, it is noted, that the quasicrystal forms at low 
temperatures, possibly reducing the amount of entropic 
disorder. We speculate, that an (aperiodic) i-TiZrNi qua- 
sicrystal, or a very similar large-cell crystal, could be a 
ground state phase. This idea is investigated in a subse- 
quent paper. 



V. CONCLUSION 

We have developed the first realistic structural model 
for the icosahedral Ti-Zr-Ni quasicrystal. Our model is 
based on atomic decorations of the canonical cell tilings. 
For the structural refinement a new method of con- 
strained least squares analysis is used, combining X-ray 
and neutron diffraction data with results from ab initio 
relaxations. The refined structural model exhibits only 
small amounts of chemical disorder between Ti and Zr. 

This structural model should be helpful for a better 
understanding and interpretation of experimental data 
such as, for example, EXAFS data, pcT-isotherms, rela- 
tions to crystalline phases and the ternary phase diagram. 
Furthermore this model will enable us to investigate the 
possibility of i-TiZrNi as a ground state quasicrystal. 
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APPENDIX A: DIFFRACTION FROM 
APPROXIMANTS 

In this Appendix a method is presented for the cal- 
culation of the structure factors of the icosahedral qua- 
sicrystal from large cubic approximants. 

Assuming that the local structure of the approximant 
closely resembles the structure of the quasicrystals, the 
quasicrystal structure factors can be approximated by the 
structure factors of a large periodic approximant. How- 
ever, a direct comparison of the structure factors of the 
approximant to measured structure factors of the icosa- 
hedral phase is not possible due to the different symme- 
tries of the structures. This leads to different multiplic- 
ities of corresponding reflections for both structures and 
a splitting of orbits of icosahedral reflections in the case 
of the approximant. The problem is further complicated 
by the fact that the Bragg reflections of the quasicrystal 
form a dense set, which is not the case for the periodic ap- 
proximants. Inequivalent low intensity reflections, there- 
fore, can become degenerate in the approximant. The 
problem of degenerate reflections in the approximant be- 
comes increasingly important in smaller approximants. 
Since, for structural investigations only Bragg reflections 
with intensities above a threshold are used, the degen- 
eracy problem can be avoided by using a large enough 
approximant. 

The six dimensional representation of icosahedral qua- 
sicrystals enables one to relate their diffraction patterns 
to the diffraction patterns of their related crystalline ap- 
proximants. The reciprocal lattice (RL) of the approxi- 
mant is obtained when the golden mean, r, which appears 
in the icosahedral basis of the physical space, is replaced 
by the rational number r„ = f n +i/fn- The mapping of 
the icosahedral RL to the cubic RL is induced by keep- 
ing the same six indices but using the rationally related 
basis vectors. Every icosahedral reciprocal lattice vector 
(RLV) maps to a unique cubic RLV. The maximum cu- 
bic subgroup of the icosahedral symmetry group, m35, is 
given by m3. While the order of the icosahedral group 
is 120, the order of the cubic subgroup is only 24. Thus, 
the image of one orbit of equivalent icosahedral reflec- 
tions is the union of up to five cubic orbits. Although 
an icosahedral RLV vector does not map to a cubic RLV 
in the same direction, any icosahedral RLV in a given 



fundamental domain of the cubic group maps to a cu- 
bic RLV in the same fundamental domain. Thus, to get 
the standard representative cubic RLV's corresponding 
to one icosahedral RLV orbit, we only need to check the 
images of all icosahedral RLV's in that orbit whose orien- 
tations fall in the standard principal triangle of the cubic 
group. 

The calculation of the structure factors for the icosa- 
hedral phase from the cubic approximant proceeds in two 
steps. First, for every orbit of icosahedral reflections with 
complementary components smaller than some cutoff, 
Q^ ax , the multiplicity, M; CO) is determined by analyzing 
if the representing RLV points into any high-symmetry 
direction of the icosahedral point group. The icosahedral 
RLV's of the orbit which fall into the fundamental do- 
main of the cubic point group are determined. Depend- 
ing on the multiplicity of the icosahedral orbit this results 
in a set of up to five RLV's. For each of these icosahedral 
RLV, the corresponding cubic RLV in the approximant is 
determined by the projection method, where the golden 
mean, r, in the basis vectors of the physical space is re- 
placed by the rational approximant, T n . The multiplicity, 
M a , of the orbit represented by each of the cubic RLV's 
is calculated by analyzing if the cubic RLV points into 
any high-symmetry direction of the cubic point group. 

In the second step, to approximate the structure fac- 
tor of the icosahedral RLV, the structure factors of the 
corresponding cubic RLV's of the periodic approximant 
are calculated and averaged. This average could be taken 
over either the amplitudes or the intensities of the struc- 
ture factors. While for reflections with a high intensity, 
both methods of averaging result in nearly the same val- 
ues, for low intensity reflections the right choice of av- 
eraging is more important. Since all the cubic approx- 
imants are centrosymmetric, the structures can be cen- 
tered such that the phases of the structure factors are 
either or tt, resulting in a plus or minus sign for the 
structure factors. The correct centering of the approx- 
imant structures involves a shift in the six-dimensional 
space. The 1/1 approximant structure requires no shift. 
For the cubic 2/1, 3/2, 8/5 and 13/8 approximants of the 
canonical cell tiling, the six-dimensional centering vector 
is given by 4(1,0,1,1,0,0). These approximants have 
the symmetry group Pa3. That means that there is one 
3-fold axis which passes through the node at the origin in 
real space. That 3-fold axis is chosen to be in the cubic 
(111) direction which corresponds to the six-dimensional 
vector (1,0,1,1,0,0). 

For reflections with nearly vanishing intensity, the 
splitting of the reflection in the approximant leads to 
a distribution of structure factors around zero. An av- 
eraging over intensities would lead to a positive value, 
while for averaging over amplitudes the resulting value is 
closer to zero and provides a better approximation. In 
this work, the icosahedral structure factor is estimated by 
averaging over the amplitudes of the corresponding struc- 
ture factors of a large cubic approximant (see Eqs. (^) 
and (pi)). In order for the approximated structure factors 
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to closely resemble the icosahedral structure factors, the 
atomic structure of the approximant has to resemble the 
structure of the icosahedral phase closely, which is more 
likely the case for larger approximants. 

APPENDIX B: AVERAGING OVER 
AMPLITUDES OR INTENSITIES 

To compare the accuracy of the prediction of the struc- 
ture factors from averaging over amplitudes and over in- 
tensities the method is applied to the Penrose tiling, since 
for the ideal quasiperiodic Penrose tiling, the structure 
factors can be easily calculated by Fourier transformation 
of the acceptance domain. If the vertices of the Penrose 
tiling are identical scatterers, the structure factors of the 
ideal quasiperiodic Penrose tiling are determined by the 



Fourier transformation of the triacontahedral acceptance 
domain l3 

In Figure || the structure factors for the ideal quasiperi- 
odic Penrose tiling are compared to the estimated val- 
ues from the cubic Penrose 1/1 and 3/2 approximant 
structures. While for the 1/1 approximant the struc- 
ture factors deviate strongly at lower intensities, already 
for the 3/2 approximants the estimated structure factors 
closely resemble the icosahedral structure factors over 
about three orders of magnitude. Averaging over ampli- 
tudes instead of intensities leads to significant improve- 
ments in the case of the 1/1 approximant, but only minor 
changes for the 3/2 approximant. It is noted here, that 
the Penrose tiling decorated with unit scatterers on the 
vertices is likely to oversimplify the problem and that in 
more realistically decorated structures, larger deviations 
could occur. 
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FIG. 1: Structure of the Bergman cluster consisting of a cen- 
tral Ni atom (left), a small icosahedron of Ti atoms (center) 
and a large icosahedron of Ni atoms with Zr atoms on the 
face centers (right). 



FIG. 2: The four canonical cells. The b-bonds and c-bonds 
are represented by the double and single lines respectively. 
The face names are also indicated. 



FIG. 3: Decomposition of the canonical cells tiling into pro- 
late and oblate rhombohedra as well as the rhombic dodeca- 
hedron. 



FIG. 4: Decoration of (a) the rhombic dodecahedron, (b) the 
prolate rhombohedron and (c) the oblate rhombohedron. The 
right-hand pictures show the same decorations with interiors 
visible. 



FIG. 5: Decoration of (a) the PR Y (C), (b) the PR D , (c) the 
PRy(D) and (c) the OR. 



FIG. 6: Comparison of the experimental X-ray diffraction 
and neutron scattering data to the simulated diffraction of 
the icosahedral structure. 



FIG. 7: Partial pair distribution functions of the W-phase 
(top) and the 8/5-approximant model (bottom). 



FIG. 8: Comparison of the structure factors of the icosahe- 
dral Penrose tiling to the estimated values from the periodic 
1/1 and 3/2 approximants. The results from averaging over 
intensities is show in the top panels and for averaging over 
amplitudes in the bottom panels. 
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